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Abstract 



A nonlinear generalization of the Fluctuation-Dissipation Theorem (FDT) 
for the n-point Green functions and the amputated 1PI vertex functions at 
finite temperature is derived in the framework of the Closed Time Path for- 
malism. We verify that this generalized FDT coincides with known results 
for n = 2 and 3. New explicit relations among the 4-point nonlinear response 
and correlation (fluctuation) functions are presented. 
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I. INTRODUCTION 



Within the framework of non-equilibrium statistical mechanics a successful linear re- 
sponse theory has been established by Kubo |l| and other authors [@,[|. One of the focuses 
in linear response theory is the Fluctuation-Dissipation Theorem (FDT) |2| which describes 
the relation between the response and the correlation (fluctuation) functions. To seek for a 
nonlinear generalization of the FDT is an important and interesting issue [PHTT|j in nonlinear 
response theory which deals with functions of many time arguments. Earlier attempts 
to derive a nonlinear generalization of the FDT were based on studies of the KMS condition 
fUU, the condit ions of time reversal invariance, and the exact algebra of transformations 
between different representations of the thermal Green functions (see Eq. (|2l| ) below). As 



stated in Ref. ||10|| , these attempts failed since none of these relations by themselves provide a 
nonlinear generalization of the FDT, and their correct combination has not yet been found. 

It was shown in Refs. flJTOl that the Closed Time Path (CTP) formalism suggested 
by Schwinger |13|] and further elaborated by Keldysh Jl4| provides a suitable theoretical 
framework for nonlinear response theory near thermal equilibrium. By combining the three 
types of relations mentioned above, the authors of Ref. |H]] arrived at a plausible nonlinear 
generalization of the FDT for the nonlinear response functions, without, however, giving an 
explicit proof of its general validity. Their suggested nonlinear generalization of the FDT 
(Eq. (5.60) in Ref. Jl0|) involves a parameter e = ±1 which depends on the signature under 
time reversal of the physical quantity considered. For arbitrary field operators it is not 
clear which value of e should be taken, which makes it difficult to verify the correctness 



of the generalized FDT given in ||10|| . In this paper we will give a complete derivation of 
the generalized FDT for nonlinear response functions. Instead of time reversal we exploit 
the closely related "tilde conjugation" transformation. We find that the correct generalized 
FDT equations contain no sign parameter e. 

Several recent investigations jl^[15]-|24[] of n-point functions (n > 3) and their spectral 
representations have tried to clarify the analytical continuation from the imaginary time 



2 



formalism (ITF) p5| to the real time formalism (RTF) |T3||14]j23|]26|| , the RTF Feyman rules 
and their transformation between various RTF representations, and the calculation of the 
temperature dependent running coupling constant in QCD. The generalized FDT plays a 
crucial role in the determination of spectral representations for the retarded and advanced 
ra-point Green functions in the CTP formalism [|TI|,^,^T|,^] . In Ref . JTI] the corresponding 



relations were derived for the retarded and advanced 3-point functions, by way of a rather 
tedious evaluation of the corresponding Lehmann spectral representation. As we will show, 
these same relations can be deduced very easily from the generalized FDT which, however, 
applies equally well to more general situations, n > 4. Due to the known intimate relation 
between transport coefficients and higher order n-point Green functions as provided by the 
Kubo formalism [P,[27|-|29|l , we expect the generalized FDT to become a useful tool for the 
calculation of collective phenomena in, say, a quark-gluon plasma. As a non-trivial example 
for its application we give explicit expressions for the relations between the 4-point nonlinear 
response and correlation functions which may be useful for a RTF calculation of the shear 
viscosity in scalar field theories ]2"7|-|2U . 



This paper is organized as follows: In Sec. [TT] we review some important relations be- 
tween various components of the rz-point RTF Green functions in the CTP formalism and 
the general expressions for the nonlinear response functions which will be used in the fol- 



lowing sections. In Sec. |T| we derive the generalized FDT for nonlinear response functions. 
By solving the generalized FDT equations, explicit relations among the correlation (fluctu- 
ation) and response functions are obtained for the cases of n = 2, 3, 4 and checked against 



known results. In Sec. [TV| we repeat the derivation for the amputated 1PI vertex func- 
tions. A short summary is given in Sec. [V[ In the appendix we deduce relations between 
the vertex functions in the physical representation and the i?/A-functions by Aurenche and 
Becherrawy ||20|| . 



II. NONLINEAR RESPONSE THEORY IN THE REAL-TIME FORMALISM 
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A. n-point Green functions and their generating functional 



In the literature there exist three equivalent representations of the RFT n-point Green 
functions and their generating functional |T0| , |30[| : the "closed time-path" , "single-time" and 



"physical" representations. 

1. In the closed time-path (CTP) representation, the n-point Green functions 
G p (l, . . . , n) (which may include disconnected parts, especially if the fields develop a vacuum 
condensate) and their generating functional Z[J(x)] are defined as 



c f (i ») = (-.)'-'W(i)'-« 

(i) 

Z[J(x)\ = (T p [exp(i J d 4 x J(x)4>(x))]). (2) 



8»Z[J(x)\ 
1 ' 5J(l)---5J(n) 



Throughout this paper the field operator <p will denote a bosonic (elementary or composite) 
field; the obvious extension to fermions is left to the reader. J(x) is a classical external 
source, and (. . .) the thermal expectation value. The numbers 1, . . . , n stand for the coordi- 
nate arguments x\, . . . , x n in 4-dimensional Minkowski space, and p denotes the integration 
contour (the closed time-path) which consists of a branch G\ running from negative infinity 
to positive infinity and another branch Ci running back from positive infinity to negative 
infinity. T v represents the time ordering operator along this time-path; it is the usual time 
ordering operator T on the branch C\ and the anti-chronological time ordering operator T 
on the branch C2. 

2. Instead of the closed time-path which covers the real time axis twice one can use a 
single real time parameter if one distinguishes the external source J\ on the branch C\ from 
J2 on the branch C2, setting them equal (and to zero) only at the end of the calculation. 
Then Eqs. (|l]) and (|2|) can be rewritten as 

G ai ... a „(l, . . . ,n) = {-iY- l (T p [^ ai {l) ■ 



SJ ai (l)---5J an (n) 



(3) 

J 1= J 2 =0 



Z[Ji(x), J 2 {x)\ = (t p exp(i f d A x [J 1 {x)4>\{x) - J 2 {x)^ 2 {x)\ 



(4) 



Here the single time parameter varies from — oo to oo, and Oi, a 2 , ■ ■ ■ , a n — 1; 2 indicate on 
which of the two branches, C\ or C 2 , the fields are located. The definition (|j) gives, for 
example, 



^2^(1, • • • , n) = H)"- 1 (T[0(1) • ■ • 0(m)]T[0(m + 1) ■ ■ ■ 0(n)]) . 

m n — m 

3. The "physical" or r/a representation of the formalism is defined by setting 



(5) 



J a (x) = Ji(x) - J 2 (x) , J r (x) = ~(Ji(x) + J 2 {x)) , 

4> a (x) = $x(x) - 4> 2 (x) , 4> r (x) = ~(J>i{x) + 4> 2 (x) 



(6a) 
(6b) 



Eqs. (|3|) and (f|) are then transformed into 



G Ql ... a , l (l,...,n) = H^ 1 2 



"- l ' 1 "'- 1 W ai (l)..^ 0l .(n)]> 



. ir on,-l ^^[J a (x),J r (x)] 

v J 5J Sl (l)---5J fin (n) 



t/ a — 1/7- — 



(7) 



Z[J a (x), J r (x)] = {Tp expfi / d 4 x [J a (x)(f> r (x) + J r (x)0 a (x)] 



Here ai, . . . , a n = a, r, and n r is the number of indices r among (oti, a 2 , . . . , a n ). For ctj = a 
one defines <5j = r, and vice versa. One shows generally that [p~0|,pOj| 



G aa ... a (l, . . . ,n) = 0. 



(9) 



Furthermore, the functions G raa _ a , G ara ^ a , G aar ^ a , . . ., with only one index r can be ex- 
pressed |T0|j30[| as sums over expectation values of n—\ nested commutators which are just 



the fully retarded n-point Green functions defined in Ref. [ 3~I] . G> r ... r is a sum over expecta- 
tion values of n — 1 nested anticommutators; it is the n-point correlation function, with no 
retarded or advanced relations among any of its time arguments. The remaining components 
Gai...a„ (1) ■ ■ ■ ,n) are sums over expectation values of combinations of n — 1 nested commu- 
tators and anticommutators, with both retarded and advanced relations between some of 
the n time arguments. 



The generating functional for n-point connected Green functions in the closed time-path 
representation is defined as 



W[J{x)\ = -i\nZ[J(x)} . 



(10) 



In analogy to Eqs. (|J) and (^), we can define its single-time representation W[Ji(x), ^{x)] 
and its physical representation W [J a (x), J r (x)]. Replacing in Eqs. the generating 

functional Z by W we obtain the connected Green functions Gp(l, 2, . . . , n), G c ai an (1, . . . , n) 
and G c ai ... Qn (l, • • • ,n) in the three different representations. 

Finally, the amputated 1PI n-point vertex functions 1^(1, . . . ,n) and their generating 
functional r[</?(a;)] in the closed time-path representation are defined as 



rf (l,...,n) 



ir ^r[y(x)] 
' 8ip{l)---5ip{n) 



<t=(4>) 



T[(p(x)] = W[J(x)\ - / d 4 x J(x)<p(x) 



with the classical field 



(p(x) 



SW[J{x)\ 
6J(x) 



In the single-time representation these definitions are rewritten as 

1 ai...a n v ' ■■■i lb ) — \ r /n % x , 

/co , 
d 4 x ( Ji(x)(/?i(x) - J 2 (2:)y?2(V) 
-co v 

with the classical fields on the two branches of the closed time-path 

SWiJ^x), J 2 (x)] _ WfJiCx), J 2 (x)] 



y?i(x) 



5Ji(s) ' 5J2(x) 

The corresponding definitions in the physical representation are 



r (n) d n) 

1 ai...ain v- 1 -' ■■■}"') 



-l) n 2 



! ^ n r[y a (a:),y r (x)] 



ipa=0,<p r = (4>) 



r[(p a (x),tp r (x)\ = W[J a (x), J r (x)] - [ d^X (j a {x)ip r {x) + J r (x)^ a (x)) 



(11) 

(12) 



(13) 



(14) 
(15) 



(16) 



(17) 
(18) 



where n a is the number of indices a among (aj., . . . , a n ), and 

(p a (x) = ipx(x) - <p 2 (x) , yv(ar) = ~{<fi{x) + (f 2 ( 
Corresponding to Eq.(Q) we have 



x 



(19) 



rW r (l,...,n) = 0. 



(20) 



The relation between the r/a vertex functions a and the R/A vertex functions 
defined in Ref. [2(J can be found in Appendix [A]. 



In Refs. |T0| , p^j30|| the following transformation law between the single-time and physical 
representations for the n-point Green functions (including disconnected parts) was estab- 
lished: 



G ai ... an (l, ...,n) =2 2 G ai ... an (l, . . .,n)Q aiai ---Qc 



(21) 



Here 



Qal — —Qa.2 — Qrl ~ Qr2 



V2 



(22) 



are the four elements of the orthogonal Keldysh transformation for 2-point functions |L4| . In 
the following, a summation over repeated indices is understood. In a similar way we derive 



^o 1 ...o n (^' 



,n 



2 2 Gai...a„ (!)•••) n )Qaiai " ' ' Qa n a n 



itU (1, ■ • • , n) = 2 1 -«rg). ^ (1, . . . , n)Q aiai ■ ■ ■ Q anan 



(23) 
(24) 



The transformation laws for the connected n-point Green functions and the amputated 1PI 
vertex functions are thus identical to those for the normal n-point Green functions (including 



disconnected parts), contrary to a remark made on page 24 in Ref. [10 



B. Tilde conjugation and KMS condition 



The tilde conjugation operation is defined by reversing the time ordering in coordinate 
space j2(J. For example, the tilde conjugate of Eq. @ is 
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0^2 iu (1, - - - , n) = (-i)*" 1 (T[0(1) • • • 0(m)]T[0(m + 1) • • • 0(n)]> . (25) 

m n — m 

In momentum space one can show that tilde conjugation is equivalent to complex conjuga- 
tion: 

G ai ... an (ku ...,kn) = (-lr 1 ^...^, . . . , k n ) , (26a) 
G ai ... an (h, ...,kn) = (-lT-'Gl^Jh, ...,k n ). (26b) 

The factor (— l) n_1 here results from the corresponding factor (— i) n ~ l in the definition 
(0). Energy-momentum conservation requires k\ + fc 2 + ■ ■ ■ + k n = 0. The relations ( [26]) 
carry over without modification for the connected Green functions G c ai an (ki, . . . , k n ) and 
G c ... a (ki, k 2) ■ ■ ■ , k n ). For the amputated 1PI vertex functions we have 

r£Lj*i, ...,*„) = -Itilfe • • • A) , (27a) 

r£U(*i, ...,*») = -rSr^(*i, ...,k n ). (27b) 

Using the KMS condition for periodicity of the Green functions in imaginary time |],[12| one 
derives in coordinate space 

G Bl ... an (l,...,n) = exp li/3 d u \ G ai ...a n (l, ■ ■ ■ ,n) , (28a) 

\ {i\a i= 2} J 

fS.. a „(l,"--,n) = exp(^ £ aJr^U(l,...,n), (28b) 

V {i|a i= 2} / 

where /3 is the inverse temperature, and = 2,1 for a, = 1,2. The connected Green 
functions C?£ ... (1, • • • ,n) satisfy the same relation Eq. ( |28a|) as the disconnected ones. In 
momentum space the exponential prefactors become products of simple inverse Boltzmann 
factors e 13 ^ . 

C. General expression for nonlinear response functions 



Let us consider a thermal system with global equilibrium initial conditions at to = — oo 
which is driven out of equilibrium by a real time-dependent external source J(t) which 
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is switched on adiabatically. In the closed time-path representation we define the Green 
functions for a nonvanishing external source by 

D(l n) =»(-!)" 6nZ[J(x)] (29) 

As shown in Refs. P,[T0||, the observables in nonlinear response theory are the correlation 
functions D r , D rr , D rrr , ... in the physical representation. These functions can be expressed 
as a Volterra series in powers of the external source J, with G ra , D raa , D rra , ... as their 
coefficients. Taking J\{x) = ^(x) = J{x), the nonlinear response can be expressed as 

D r {l) = G r {l) - Jd2 G ro (l, 2) J(2) + ^Jd2d3 G raa {l, 2, 3) J(2) J(3) + . . . , (30a) 
D„(l, 2) = G rr {l, 2)-Jd3 G rra (l, 2, 3) J(3) 

+i|rf3rf4G rraa (l,2,3,4)J(3)J(4) + ... , (30b) 
A. rr (l,2,3) = G rrr (l, 2, 3) - J c/4G rrm (l,2,3,4)J(4) 

+ 2l / M db 2 ' 3 ' 4 ' 5) J(4) J(5) + ' • • • (30c) 

Here di is a shorthand for d A X{. Following Refs. the correlation functions D r , D rr , 

are called fluctuation functions in the non-equilibrium and 
equilibrium state, respectively. The fully retarded Green functions G ra , G raa , G raaa , etc., 
correspond to linear, second-order, third-order, and higher order response functions of the 
averaged physical field cp to the external source J. Similarly, G rra , G rraa , G rrra , . . . are the 
response functions of the fluctuations at different order. 



III. NONLINEAR GENERALIZATION OF THE FLUCTUATION-DISSIPATION 
THEOREM FOR n-POINT GREEN FUNCTIONS 

Eqs. ( |28a| ) show that not all 2 n components of the n-point Green functions G ai ... an in the 
single-time representation are independent of each other. Taking also into account Eq. ( pip 
with Eq. (|9|), one sees that there are at most 2 n_1 — 1 independent components. We will 
soon see that the same is true for the physical representation of the n-point Green function, 
C 
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As a simple example consider the 2-point Green function. According to @ G aa = 0. 
The remaining three components are defined in coordinate space by 

G ro (l,2) = G rct (l,2) = -^ti-t 2 )([<Kl),0(2)]>, (31a) 
G ar (l, 2) = G adv (l, 2) ee i9(t 2 - h)([4>(l), 0(2)]) , (31b) 
GW(1, 2) = G cor (l, 2) ee -»<{#(1), 0(2)}) . (31c) 

In momentum space the retarded and advanced 2-point functions, G vet and G adv , are related 
by complex conjugation, while the correlation function G cor is related to them by the well- 
known fluctuation-dissipation theorem 0, in our notation 

G rr (k) = (l + 2n(k )) (G ra {k) - G ar {k)) . (32) 

Here n(k°) = [e^ fe ° — 1] _1 is the Bose-Einstein distribution function. The FDT thus plays 
an important role for the analytic structure of the thermal 2-point function. 

We will now derive similar relations among the n-point (n > 3) retarded and advanced 
Green functions in the physical representation. Since, as discussed in the preceding section, 
these functions are related to the higher order fluctuations and nonlinear response functions, 
the relations among them contain the desired nonlinear generalization of the FDT. 

Starting from the transformation (|2lD and using Q rdji = Q rai and Q aai = —Q aai we get 

(1, ...,n) ±G Sl ...a„(l, ...,n) (33) 
= 2 1 "? ^ (l ± (-l) n ^-^)G ai ... an (1, . . . , n)Q aiai ■ ■ ■ Q anan , 

ai,...,a n =a,r 

where n a (ai, . . . , a n ) counts the number of a indices among (oti, . . . , a n ). Eqs. (|33|) form 
2«-i independent equations. We now consider separately the contributions to the sum on 
the r.h.s. of (|33| ) with fixed even and odd numbers of a indices n e a and n°, and define 

Gnl,(ai...a n )0-i ■ ■ ■ > n ) = G ai _ an (l, . . . , Tl)Q aiai ■ ■ ■ Qa n a n , (34a) 

a 1 ,...,a n — a,r 
«a(o! 1 ,...,an) = n| 

(1, . . . ,n) = 22 G ai ... an {l,...,n)Q aiai ---Q anan . (34b) 

n a (ai,...,an) — n° 
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Due to Eq. (§) neither n e a nor n° a can be larger than n — 1 . Eq. (^) then gives 



Y G n%,{ ai ...a n )0-, ■ ■ ■ ,n) = 2 n 2 G ai ... an (l, . . . + G Sl ...a n (l, 



0<n^<n-l 



n) 



Y G ngi(ai ... an )(l,. . . ,n) = 2 n 2 [G ai ...a n (l,- • .,n) - G 5l ...a n (l, . . . , n) 

l<n°<n-l 

With the help of the KMS condition (|28a|) and some further algebra we find 



Y (pn°,(ai...a n )(l, ■ ■ ■ ,Tl) + <5 n e )(ai ... 0n )(l, . . . , n) 

.4 



0<n|<n-l 
= coth 



Kn»<n-1 



X! (G' ri °,(a 1 ...a n )(l, . . . , 7l) + G ? ng ,( 0l ... a „) (1, • • • , Tl) 



Y ( G <,(a 1 ...a n )(l ) • • • ,n) - G„ S ,( ai ...a n )(l, ■ ■ ■ ,U 



0<n e „<n-l 



tanh 



.4 



X! (^ng,(oi...a„)(l) •••,«)- (5 n o )(ai ... an )(l, 



Kn°<n-1 



where 



(35a) 
(35b) 



(36a) 



(36b) 



A = i(5 Y d u, 

{i|oi=2} 



(37) 



Going to momentum space and using ( [26a|) we can separate the real and imaginary parts of 
these equations. For odd values of n we obtain 



Re 



Y Gnt,,(ai...a n )(kl, ■ ■ ■ , k n ) 
0<n=<n-l 



Im 



coth ^ Y k °i Re 
Y, G n e^ ai an ^(ki, . . . , k n 



Y< Gn%,(a 1 ...a n )(ki, . . . , k r < 



0<n=<n-l 



tanh (| Y ^°] Im 

^ {i\a t =2} 



Kn°<n-2 



while for even n we find 



Im 



Y, G>i£,(ai---a„)(&l> • ■ ■ i kn) 
0<n e n <n-2 



coth[§ Y ^°] Im 

/ {*k=2} 



Yj ^ng,(ai...an)(^l) ■ ■ ■ j K, 



Y ^ng,(ai...a„)(^l) • • • j K 



Kn°<n-1 



(38a) 



(38b) 



(39a) 
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Re 



Gng,(a 1 ...a n )(k 1 , . . . , k n ) 



0<nf<n-2 



tanh 



{i\a,i=2} 



z2 G n °,( ai ...a n )(ki, . . . ,k n ) 



\<n°<n-l 



(39b) 



Eqs. ( |38| , |39| ) are the nonlinear generalization of the FDT for n-point Green functions. Each 
of these two sets of equations contains 2 n ~ 1 pairs of (real) relations among the components 
G ai ... an in the physical representation. Taking also into account Eq. (£5|) we see that at most 
2 n_1 — 1 (complex) components are independent, as promised above. 

We should point out that the nonlinear generalization of the FDT for the connected n- 
point Green functions is the same as Eqs. (|38| , |39|) . This is because G and G c have the same 

) and satisfy the same KMS condition (|28a|) . The 



transformation properties (21) and (| 
following examples will therefore be written down only for the standard n-point functions 
which may contain disconnected contributions like, e.g., vacuum condensates. 



A. 2-point function 



Combining Eqs. @ and (21) for n = 2 leads to 



G22 — ~^{G rr — G ra — G ar ) , 
G12 = -z(G rr — G ra + G ar ) , 
We can then write down the definitions 



G21 — —(G rr + G ra — G aT/ 

Gl\ — ^(Grr + G ra + G ar ) 



(40a) 
(40b) 



for the two relevant cases n e „ = 0, n° = 1: 



G 



G r 



0,(22) — ^rr , 
^0,(21) = G rr , 



£"1,(22) — —G ra ~ G ar ] 
Gl,{21) = G ra ~ G ar . 



Substitution into fl39|) yields in momentum space 



Im 



G ra (h, k 2 ) + G ar (ki, k 2 ) 



Re G rr (ki, k 2 ) = 



Im G rr (ki, k 2 ) = coth — — Im G ra (ki,k 2 ) 

Bk° r 

ReG rr (A;i, k 2 ) = tanh— -i- Re G ra (ki, k 2 ) 



G ar (ki, k 2 ) 
G ar {ki, k 2j 



(41a) 
(41b) 

(42a) 
(42b) 
(42c) 

(42d) 
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The zeros in the first two equations follow from energy conservation, k® + k® = 0, which 
leads to coth[/3(fc° + fc°) /2] -> oo and tanh[/5(fc° + fcg) /2] -> 0. Eqs. (g2|) are easily solved by 



G> r (&) = coth -— |G> a (fc) - G ar (k) 



G ra {k) = G* ar (k) 



(43) 



The first of these equations is the well-known fluctuation-dissipation theorem (]32[) . 

Substituting Eq.fl4"3|) into Eq.(40) and using the two-component column vectors intro- 



duced in Refs. p3|j26[| one can express Eq. (|40| ) in terms of an outer product of these column 



vectors 23,11 



G aia2 (k) - G ra (k) Q ( + J 0) A G ar (k) ( + ( ^ 0) 



(44) 



B. 3-point function 



For n = 3 Eqs. @ and (pTT ) give 



G222 = 


(G rrr 


G rra 


Gf-ar 


G arr 


"I - G r aa 


G ara 


G a ar) ; 


(45a) 


G211 = 


— (G rrr 


-\- G rra 


G rar 


Garr 


~t~ G r aa 


G ara 


G a ar) j 


(45b) 


G121 = 


— (G rrr 


-\- G rra 


G ra r 


~)~ G a rr 


Graa 


G a ra 


Gaar) j 


(45c) 


C112 = 


— (G rrr 


G rra 


G r ar 


G a rr 


Graa 


G a ra 


G a ar) j 


(45d) 


C221 — 


— (G rrr 


-\- G rra 


Grar 


Garr 


Graa 


Gara 


G a ar) j 


(45e) 


G212 = 


— (G rrr 


G rra 


G r ar 


G a rr 


Graa 


G a ra 


G a ar) j 


(45f) 


G\22 = 


—{G rrr 


G rra 


Grar 


G a rr 


G T aa 


Gara 


G a ar) j 


(45g) 


Gin — 


—{G rrr 


-\- G rra 


G ra r 


G a rr 


G r aa 


G ar a 


G aa r) ■ 


(45h) 



The definitions (|34"D yield 



^0,(222) — Co,(211) — Co,(121) — ^0,(112) — G rrr ', 



Cl,(222) — —Grra ~~ G rar ~ G a 
Cl,(211) = G rra + G rar ~ G arr 



G 



2,(222) 



G r aa "I - G a ra "I - G a ar i 



^2,(211) — C raa — G ara — Gaar 
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(46a) 
(46b) 
(46c) 



Gi, 
Gi, 



(121) 



(112) 



G r ra G r ar G ar r , G 

G rra -\- G rar -\- G arr , 



2,(121) 



G 



2,(112) 



~G raa -\- G ara G aar ; 
G r aa G ara -\- G aar 



(46d) 
(46e) 



Substituting these quantities into Eqs. (|38D we get in momentum space 



Re 
Im 
Re 
Im 
Re 
Im 
Re 
Im 



G r ra G r ar G arr , 
G r rr G r aa G ara -\- G aar 
G rrr ~t~ G ran. G ar a G aar 



G r rr G r aa G ara G aar 



G r rr G r aa G ar a G aar 



Grrr G r aa G ar a G a 
G r rr G r aa G ar a ~\~ G a 
G r rr G r aa G ar a G a 



0; 

coth — — ^ Re 

tanh ^— ^ Im 
2 

coth — - Re 
2 

tanh — - Im 
2 

8k° 
coth — Re 
2 

tanh im 



Grra G r ar G arT 
G r ra G r ar G ar r 
G r ra G T ar G arr 
G r ra G r ar G arr 
G r ra G r ar G a 
G r ra G r ar G c 



(47a) 
(47b) 
(47c) 

(47d) 
(47e) 
(47f) 
(47g) 
(47h) 



Again the zeros in the first two equations result from energy conservation. To solve Eqs. (fI7|), 
we choose G raa , G ara and G aar as the 2 n ~ 1 — 1 = 3 independent physical components; 



in the notation of Ref. (TT| these are the three retarded 3-point functions Vm, Tr, and 
Tr , respectively. We further introduce the shorthands [TT) Ni = N(k®) = coth(/%j°/2) = 
1 + 2n(k®) = 1 + 2rii which satisfy the identity 



iViiVa + N 2 N 3 + N 3 N ± 



1 . 



(48 



Solving Eqs. ( Ww in terms of the three selected components we find 

G rar = ^l{G* ar a ~ G aar ) + N 3 (G* ara ~ G raa ) , 
Garr = ^2(G* aa — G aar ) + N 3 (G* aa — G ara ) , 
G rra = ^l[G aar — G ara ) + ^2{G* aaT — G raa ) , 
G rrr = G* aa + G* ara + G* aaJ . + N 2 N 3 (G r aa + G* aa ) 

+N 1 N 3 (G ara + G* ara ) + N^Gaar + G* aar ) . 



(49a) 
(49b) 
(49c) 

(49d) 
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Noting the identities G rar = Tp, G arr = Tpi, G rra = Tp Q , and G rrr = Tp between the 
3-point functions m r/a notation and the 3-point functions used in ]TTJ, Eq. ( pE9"|) are seen 
to agree with Eqs. (33) in fln}| . Obviously, the algebra leading from (^) to fl49| ) is much 



simpler than the method of derivation presented in JTI] . 

Eq. ([49d|) expresses the fluctuation function G rrr in terms of the second order non-linear 



response functions (G raa , G ara , G aar ), i-e. in terms of the fully retarded vertex functions. 
Furthermore, Eqs. (f|9|a-c) express the linear response functions for the fluctuation (G rra , 
Grar, G arr ) in terms of the same set of retarded Green functions. These four relations form 
the nonlinear generalization of the FDT for the case n = 3. 

When substituting Eq. ( f49|) into (0), Eq. ([45]) can be written in terms of outer products 
of column vectors as [Tl[ 



G 

~^~G ara 

~\~G nar 



raa {l) (l + n 2 ) (l +n 3 , 

i+njyijyi+mj 

l + m \l + n 2 \l 



n 



" aar 



(1 + n 2 )(l + n 3 ) - n 2 n 3 
1 + ni)(l + n 3 ) - nin 3 
1 + rii)(l + n 2 ) - ^in 2 



(l + nj(l) (l, 
JllJll + rJ 



(50) 



C. 4-point function 



Following the same steps as for n — 2, 3, we get from Eqs. (|J), fl2"T|), and (|4j) 



^0,(2222) — Co,(2221) — ^0,(2212) — ^0,(2122) — ^0,(1222) 



G 



0,(2211) — Go,(2121) 



Gq. 



(1221) — G rrrr ; 



(2222) 



G 



~G rrra G rrar G rarr G arrr , 



2,(2222) G rraa -\- G rara -\- G raar -\- G arra -\- G arar -\- G aarr , 



G 



3,(2222) 



~G raaa G araa G aara G aaar , 



(51a) 



(51b) 



*~ T 1,(2221) ^rrra ^rrar ^rarr ^arrr j 
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G 2 
G 3 
G l 
G 2 
G 3 
Gi 
G 2 
G 3 
G 1 
G 2 
G 3 
G x 
G 2 
G 3 
Gi 
G 2 
G 3 
Gi 
G 2 
G 3 



,(2221) 
,(2221) 
,(2212) 
,(2212) 
,(2212) 
,(2122) 
,(2122) 
,(2122) 
,(1222) 
,(1222) 
,(1222) 
,(2211) 
,(2211) 
,(2211) 
,(2121) 
,(2121) 
,(2121) 
,(1221) 
,(1221) 
,(1221) 



G rraa G rara -\- G raar G arra -\- G arar -\- G aarr 

G ra aa G araa ~\~ G aara G aaar , 

G rrra -\- G rrar G rarr G arrr , 

G rraa -\- G rara G raar -\- G arra G arar -\- G aarr 

G ra aa G ara a G a ara G a aar j 

G rrra G rra r ~\~ G ra rr G arrr , 

G r raa G ra ra G r aar G a rra G a rar G aa rr ) 



Graaa G ara a ~l~ G a ara ~t~ G a aar j 



Grrra G rra r G ra rr G arr r j 
G r raa ~\~ G ra ra G ra ar G ar ra G a rar G aa rr > 

Graaa G a raa G a ara G aa ar j 
Grrra G rrar G rarr G arr r > 

Grraa G rara G r aar G arra G ara r ~\~ G aa rr i 

Graaa G a raa ~\~ G a ara G aa ar j 
Grrra G rrar ~f~ G rarr G arrr , 

G rra a G ra ra G T aar G a rra G a rar G aa rr 

Graaa G a raa G a ara G aa ar > 
Grrra G r rar G r arr G ar rr > 

G rr aa G rara ~\~ G T aar G a rra G a rar G aa rr 
Graaa G ar aa G a ara G aa ar ■ 



Substituting these quantities into Eq. (|39p, we get in momentum space 



Im[Gi, 



(2222) + (3^(2222) 



o, 



Re [Gq, (2222) + G 2 



,(2222)J 



(3k c 4 



Im [Go, (2221) + C 2 ,(2221)] = — COth — ^— Im[G f i ) (2221) + C 3) (2221)] , 
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Re[G , 


(2221) 


+ G 2 , 


(2221)] 


lm[G , 


(2212) 


+ G 2i 


(2212)] 


Re[G , 


(2212) 


+ G 2 , 


(2212)] 


Im[G , 


(2122) 


+ G 2 , 


(2122)] 


Re [G , 


(2122) 


+ G 2 , 


(2122)] 


Im[G , 


(1222) 


+ G 2 , 


(1222)] 


Re [G , 


(1222) 


+ G 2 , 


(1222)] 


Im[G , 


(2211) 


+ G 2 , 


(2211)] 


Re [G , 


(2211) 


+ G 2 , 


(2211)] 


Im[G , 


(2121) 


+ G 2 , 


(2121)] 


Re [G , 


(2121) 


+ G 2 , 


(2121)] 


Im[G , 


(1221) 


+ G 2 , 


(1221)] 


Re [G , 


(1221) 


+ G 2 , 


(1221)] 



(3k° 

- tanh Re[Gi i (222i) + G 3i ( 222 i) 

- COth Im[Gi i(2 212) + G 3) (2212) 

(3k° 

- tanh — - Re[G f lj (22i2) + ^3,(2212) 

- COth —± Im[Gi i( 2l22) + G 3) (2122) 

(3k° 

- tanh Re[Gi i(2 i22) + (^3,(2122) 

(3k° 

- COth Im[Gi i( i222) + G 3; (i222) 

(3k° 

- tanh — - Re[Gi i( i222) + (^3,(1222) 



: coth 
tanh 
coth 
tanh 

: coth 



2 






kl) 


2 






kl) 


2 






kl) 


2 




m + 


kl) 


2 




+ 


kl) 



Im[Gi i (22ii) + G 3; 
Re[Gi ; (22ii) + G 3 . 
Im[G f i j (2i2i) + G 3; 
Re[Gi ; (2i2i) + G 3i 
Im[Gi 5 (i 22 i) + G 3; 
Re[Gi ; (i22i) + G 3 . 





(52d) 




(52e) 




(52f) 




(52g) 




(52h) 




(52i) 




(52j) 


2211)] , 


(52k) 


(2211)] ; 


(521) 


2121)] , 


(52m) 


(2121)] ; 


(52n) 


1221)] , 


(52o) 


(1221)] • 


(52p) 



Once more the zeros in the first two equations follow from energy conservation. For n = 4 
there are at most 2 n ~ 1 — 1 = 7 independent components of G aiot2a30t4 ] we choose them as 
G r aaai G araa , G aara , G aaar , G arra , G arar and G aarr . Introducing the notation 

Ni + Nj _ 1 + m + rij 



N, 



(ij) 



(53) 



< fcI > N k + N t l + n k + m' 
we can express the solution of the above equations as 

G m = ~N 2 N 3 N A G raaa + (N 2 N 3 N, + N 2 + N 3 + N 4 )G* raaa + N^N^araa 

+N 2 (N$§ + N 2 4 N$)G* araa + iVW 2 iV 4 G aara + N 3 (n$ + N 2 2 N$)G* aara 
+N 1 N 2 N 3 G aaar + N A (N$ + N 2 N^)G* aaar + A^G^ + N 2 N 3 N$>G* a 
+N 1 N 3 G arar + N 2 N 4 N$G* arar + N^G^ + N 3 N 4 N$G* aarr , (54a) 
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* arra 



Grrra ~ N 2 N 3 G raaa ~ N 2 N 4 N^f ) G* araa - N 3 N 4 N^G* aara - (N(2% + N 3 N (34)) ^loar 

-NiG arra - N 2 N ( ) G* arar - N 3 N^G* aarr , (54b) 
G rrar = N 2 N 4 G raaa - N 2 N 3 N$G* araa - (JV§§ + NlN^)G* aara - N 3 N 4 N$G* aaar 



—N 2 N^G* arra — NiGarar ~ N 4 N^ G* aarr , 



G rarr = N 3 N 4 G raaa - (ivgg + N 2 4 N$)G* araa - N^N^G^ - N 2 N 4 N$G* 



-N 3 N^G* arra - N 4 N^G* a 



'(23) 



(13), 



(24)' 



Gf«rr = (l + ^2 ^ 3 + AW + N 3 N 4 )G* raaa - N 3 N 4 G araa - N 2 N 4 G aara - N 2 N 3 G a 



n 

^raar 
G r ara 
G rrn.a 



N 4 G arra N 3 G araT N 2 G a arr 

—N 4 G raaa + N 3 N^G* araa + N 2 N^G* aara - NiG aaar + N^G* arra , 
-N 3 G raaa + N 4 N^G* araa - NiG aara + N 2 N^G* aaar + N^G* ar 



-N 2 G raaa - NiG araa + N 4 N,^G* + N 3 N,3\G* + N ( 3,G* a 



(54c) 

r 

(54d) 

r 

(54e) 
(54f) 
(54g) 
(54h) 



(24) aaar ' (24) arar ' 

r(12) (-1* . t,t(12)^y* 
'(34) ^ aara 2, 3 lV (34) °aaar "r iv (34) "aarr ■ 

This solution expresses the fluctuation function G rrrr , the linear response functions of the 
3 rd order fluctuations G rrra , G rrar , G rarr , G arrr , and 3 of the non-linear response functions of 
the 2 nd order fluctuations (G raar , G rara , G rraa ) in terms of the fully retarded 4-point vertex 
functions (3 rd order non-linear response functions G raaa , G araa , G aara , and G aaar ) and the 
remaining 2 nd order response functions of the 2 nd order fluctuations (G arra , G arar , G aarr ). 
Although these equations still have the structure of a fluctuation-dissipation theorem, they 
are clearly much more involved than for the simpler cases n = 2, 3 since fluctuation and 
response functions of different order simultaneously. 

Similar to Eq. (|50D we can express the 4-point Green function in the single-time rep- 
resentation in terms of tensors constructed from outer products of 2-component column 
vectors: 



a\ 0,20,30,4 



~G r 



1 



n 2 



n 3 



^raaa 



1J \1 + n 2 J \1 + n 3 J \1 + n 4 J 



1 + n 2 )(l + n 3 )(l + n 4 ) - n 2 n 3 n 4 



n 3 \ 1 1 + n 4 
l + n 3 )\ n 4 . 



1+nJ \1 / 

■C:r)( 



n 4 \ 
l + n 4 ) 



2 G aara 
2 aara 

2 aaar 

2 aaar 
T 2 '-■'arra 

~r 2 arar 
~r 2 aarr 



l + ni r 
n 2 V 1 / 

[l + nj 
1 + wi A 

( ^ ) 

u + V 

1 + 7li A' 
n 4 \1 

( ni ) 

( m \ 
1 + nJ 



n 2 
,1 + n 2 , 

n 2 \ 

l + n 2 ) 

" l + n 2 

n 4 

i + n 2 ; 

1 + ™2 , 



1 + ^3 



1 + n 4 



//:" 



C 



n 4 \i / 
\ AW 

{ ^ ) 

+ n 3 s 
n 3 , 



+ 



+ 



+ 



n-. 



t) (l) (l + n 4 , 
1 + nJUU 



l 



l + n 4 )V 
'l + n 2 
\ « 2 
' n\ \ 
.(1 + n 4 )V ^2 
l + n 2 \/ n 3 \1A N 

n 2 \ 
l+n 2 ) 

n 2 \ 
1+nJ 



n 3 V(l + n 3 ) 
1 +n 4 / 



Ho 



'^ + n 2 )\ 



n 3 

1 W 



+ 



"i" 2 «mi (23) 



A7(13) 
'// ar JV (24) 



iff 

Ir* /v (12) 

2 ^ T aarr iv (34) 



A/ n 4 \ 

n 3 \ ( n A \ 
l+n 3 [l + nj ' 



In deriving this expression we have made frequent use of the following relations which 
in global thermal equilibrium: 



(55) 
hold 



n(—x) — — (1 + n(x)j , (56a) 
n{x{)n{x 2 )n{x 3 ) + ^ + ^ + _ (56b) 

They are a special case of a general relation proved in the appendix of Ref . P^j . 



IV. GENERALIZED FLUCTUATION-DISSIPATION THEOREM FOR 
AMPUTATED 1PI n-POINT VERTEX FUNCTIONS 



One can derive in a similar way a nonlinear generalization of the FDT for the amputated 



1PI n-point vertex functions. In analogy to (|34|) we define 



r i°,(oi-..o„)( 1 ' • • • ' n ) 



^ 2 ^a^.On (1) • • • > n )Qaiai ' ' " Qa n a n 



aj an=a,r 

no(«l,-,"ii)=»5 



^ 2 | _ an (1, . . . , Ti)<5aiai ' ' ' Qa n a n 



(57a) 



(57b) 



aj an= a . r 

n a (ai,...,o„)=n2 



Due to ( [ZD] ) must be > 2 and n° a > 1. Using Eqs. (^) and (^) we then find instead of 



Eq. 
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V r (n) ( i 

L nS,(oi...o n )V- L ' 



V r (n) n 

1 n°,(a 1 ...a n )V- L ! • • • 



m-2 



in-2 



r (n) fi + r^ n) - (1 

1 ai...On\ i " " " ' "v ' tti...o n \- L 



r£L„(i,...,n)-r£U(i,...,n) 



Kn2<n 



With the help of the KMS condition (|28b|) and the tilde conjugation relation 
find 



Im 



2<ng<n 



/3 



coth 5 E fcHlm 



{i|a;=2} 



r y,(ai...a„)(^l' • • • > ^n) 



Kra2<ra 



Re 



H r nl,(ai...On)(^'---'^ n ) 



2<nS<n 



/3 



tanh | ^ ^ fc? ] Re 

2 {i|<H=2} 



(58a) 
(58b) 

we finally 



(59a) 



(59b) 



Yl r <,(ai...a„)( fc l' • • • ' fc « 
l<ng<n 

These relations are the nonlinear generalization of the FDT for the amputated 1PI n-point 
vertex function. Contrary to the case for the n-point Green functions, we here obtain 
formally the same relations for even and odd values of n. Eqs. (|59"D provide 2 n_1 pairs of 
(real) relations among the (complex) components of . Together with Eq. ( p0[ ) this 

leaves at most 2 n_1 — 1 independent complex components for the n-point vertex. 



A. 2-point self-energy 



Making use of Eqs. ( p0| ) and (|2l| ) in the definition (0), we find 



F (2) _ _-n(2) _ -p(2) 
1,(22) " L ra L ar 

F (2) _ F (2) _ F (2) 

L 1,(21) ra L ar i 



<2) 



r 



2,(22) 1 aa 



(2) . 



,(2) 
2,(21) 



_p(2) _ 



Substitution into (|59D yields 



Reri 2 a )(fc 1 ,A; 2 ) = 0; 

-hnVf}{hM) = coth^Im[rg(A ;i) fe) 



(60a) 
(60b) 



(61a) 
(61b) 
(61c) 
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-RelfM, k 2 ) = tanh ^ Re [rgfa, k 2 ) - V^{K K 



(61d) 



Replacing in these equations r^ 2 ^ by the 2-point Green function G and interchanging the 
indices r and a we recover Eq.(^). The solution of (^1|) is thus straightforward: 



f3k° 



r« (AO =coth^[rW (*)-!£>(*) 



r%(k) = r ( ?>(k). 



(62) 



This coincides with the results obtained in Refs. [|TI],[Tl]]. Here T$ and are the retarded 
and advanced self-energies, E ret and E adv , respectively. 

Using the column vector notation one can write the 2-point vertex function in the single- 
time representation as 



(63) 



This expression exemplifies the general substitution rule: in order to obtain the vertex 
functions T from the corresponding Green functions G, one interchanges all r and a indices 
and changes the sign in the lower components of all the column vectors. 



B. Amputated 1PI 3-point vertex 



Substituting the quantities 



.(3) 



r ( 3 ) - r 

3,(222) ~~ 1 3,(211) ~~ 1 3,(121) ~~ 1 3,(112) 



(3) 



.(3) 



_r< 3 ) 



r (3) _ r (3) _ r (3) _ r (3) p(3) r ( 3 ) , r (3) , r (3) . 

1,(222) rra 1 rar 1 arr ) 1 2,(222) raa > 1 ara ~ 1 aar > 

p(3) = r (3) , r (3) _ r (3) p(3) = p(3) _ p(3) _ p(3) . 

1,(211) rra ^ rar arr J 2,(211) raa ara aar > 

p(3) = r (3) _ p(3) , r (3) p(3) = _ r (3) , p(3) _ p(3) . 

1,(121) rra rar ^ arr J 2,(121) raa ^ ara aar ) 

r (3) = _ r (3) , r (3) , r (3) p(3) = _ r ( 3 ) _ r (3) , r (3) 

1,(112) rra ^ rar ^ arr ? 2,(112) raa ara ^ aar 



into Eq. (|59|) we find 



Im 

Re 



r(3) _j_ p(3) _j_ F (3) , r (3) 
rra ~ rar ' arr ~ aaa 



F(3) , p(3) , F (3) 
raa ' ara ' aar 



(64a) 
(64b) 
(64c) 
(64d) 
(64e) 



(65a) 
(65b) 
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Im 
Re 
Im 
Re 
Im 
Re 



r(3) _r(3) _ r (3) 

raa ara aai 

r(3) _p(3) _p(3) 
raa ara aa7 

_f(3) _i_r(3) _r(3) 

raa ' ara aai 

_r(3) i p(3) _ r (3) 
raa ' ara aai 

_F(3) _p(3) , p(3) 
raa ara ' aai 

_p(3) _ r(3) + p(3) 

raa ara ^ »™ 



coth ^-Im 



tanh 



2 



Re 



coth ^-^Im 
2 

tanh ^— ^-Re 
2 

coth— -5- 
2 

tanh — 



Im 



Re 



F(3) , p(3) _ p(3) _ p(3) 
rra ' rar arr aao 

p(3) , p(3) _ p(3) _ p(3) 
rra ~ rar arr aac 

p(3) _ p(3) , p(3) _ p(3) 
rra rar ~ arr aao 

p(3) _ p(3) _j_ p(3) _ p(3) 
rra rar ' arr aac 

_p(3) i p(3) , p(3) _ p(3) 
rra ^ rar ^ arr aac 

_p(3) , p(3) , p(3) _ p(3) 

rra rar rrrr nv7y 



arr aaa 



Replacing by the 3-point Green function G and interchanging the indices r 
recover Eq. (|47D . Choosing as independent components the fully retarded vertices 
and r^ a , we thus obtain directly the solution 

r( 3 ) = jvi fr (3) * - r (3) ) + iwr (3) * - r (3) ) 

ara H rar rra J 1 M\ rar arr J ' 

r( 3 ) = jv 9 fr (3) * - r (3) ) + iwr (3) * - r (3) ) 

raa A\ arr rra) 1 o\ arr rar J ' 

r( 3 ) = jvi fr (3) * - r (3) ) + iVofr (3) * - r (3) ) 

aar H rra rar/ 1 ^ \ rra arr/ ' 

p(3) _ p(3)* I p(3)* _i_ p(3)* , /V„A/ ('p(3) _|_ p(3)*"\ 
aaa arr ' rar ' rra ' -"J-"3V A arr ~ arr / 

+JViJV s (rS. + rfj*) + iVxiv 2 (ra + r$f) . 



(65c) 
(65d) 
(65e) 
(65f) 
(65g) 
(65h) 

and a we 
p(3) p(3) 

arri rar 

(66a) 
(66b) 
(66c) 

(66d) 



These relations are the nonlinear generalization of the FDT for the amputated 1PI 3-point 
vertex function. 

In terms of outer products of column vectors the 3-point vertex function in the single-time 
representation can be expressed as 

' 1 \ / 



1 \ / 1 



+ r ~( _ (1+ni) )(_. 



_ p(3)* 



^1 + ni)(l + n 3 ) - nira 3 



n 3 \ 

-l)(-(l + n 2 ))(-l 
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+ F (3) 



r( 3 r K(i + m)(i + n 2 ) -mn 2 



" 2 1 

i + n 2 ); v-v 

'-XV 



n 3 

-(l + *s), 



in accordance with the above substitution rule (cf. Eq. fl5D|)). 



C. Amputated 1PI 4-point vertex 



(67) 



Repeating the same technical steps for the 4-point function we find 



r (4) 



F(4) 
aara 

F(4) 
araa 

F(4) 
raaa 

r(4) 

arra 

r(4) 

arar 
r (4) 



-N 2 N 3 N 4 T<£ rr + (n 2 N 3 N 4 + N 2 + N 3 + iV 4 )rW* r + N^N^jSn- 

+N 1 N 2 N 3 T^ ra + N 4 (N$ + NiN$)r«> + N.N.T^ + N 2 N 3 N^T^> 

(68a) 



" rJv l JV 2 iv 3- L rrra ' JV 4\^ JV (24) ' Jv 3 Jv (34) J L rrra ' JV 1 JV 4 J - ra ar 

+N 1 N 3 r$ r a + N 2 N 4 N$rW; a + n.n^ + N 3 N 4 N^vi 4 > 



rraa ) 



1 aaor — iV 2 iv 3 i arrr J V 2^ 4 i V ^ 1 rarr iV 3 iv 4^ (34) J- rrar \ Iy (24) + iV 3 iV (34) J 1 rrra 



^l^raar -^2-^(24) ^roro ^3-^(34) > 

- (N[% + NlN$y&. - N 2 N 3 N$Tl?: r - N 2 N 4 N$4i 

(68d) 



)* 

rrra 



rrra 



(68b) 



(68c) 



-iv,iv; (14 ?r (4) * - Na n^t^* -^r (4) 

Jy 3 JV (23) 1 raar 1 y 4 J v (24) L rara Jv l i rrau/ 

)r< 4 >* - jV3iv 4 r (4) -A^ 2 A^ 4 r (4) -iv 2 iv",r (4) 

/ arrr o 4 rarr ^ t rrar * o rrra 



(l + N 2 N 3 + N 2 N 4 + N 3 N 4 )r^; r 
_iv 4 r (4) -^r (4) -iv 2 r (4) 

* raar o rara ^ rraa 

M r( 4 ) _i_ M aK 14 )f( 4 )* _i_ m /v( 14 M 4 )* M r( 4 ) _i_ /v( 14 )r( 4 )* 

JV4I arrr -|- l\ 3 l\( 23 j 1 rarr -|- iV 2 -'V(23) 1 rrar JVli rrra ' Jv (23) 1 raar 5 

(13) r (4)* _ AT p(4) , at n/v (13) r (4)* , /v (13) r (4)* 

(24) 1 rarr 1 v 1 1 rrar ~ J v 2 1 v (24) 1 rrra ' 1 v (24) 1 rare 



iv 2 i arrr ivii mrT -|- iV4iV( 34 j 1 rr , ar 



/v (12) r ( 4 )* ^ /v (12) r ( 4 )* /v (12) r ( 4 )* 
JV (34) 1 rrar + JV 3 JV (34) 1 rrra ' JV (34) 1 rraa ■ 



(68e) 
(68f) 
(68g) 
(68h) 



These relations are the nonlinear generalization of the FDT for the amputated 1PI 4-point 
vertex function. They can again be obtained by substituting in Q54]) G by r^ 4 ^ and inter- 
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changing r with a. The corresponding column vector representation is obtained by applying 
the above-mentioned substitution rule to Eq. (|55"|). 



V. CONCLUSIONS 

Using the transformation between the single-time and physical representations of real- 
time thermal n-point functions, the tilde conjugation relation and the KMS condition, we 
have derived the nonlinear generalization of the fluctuation-dissipation theorem for real- 
time n-point Green functions and 1PI amputated vertex functions at finite temperature. 
This generalized FDT is formally identical for disconnected and connected n-point Green 
functions. The FDT for the amputated 1PI vertex functions is obtained from that for the 
Green functions by interchanging the indices r and a. 

The generalized FDT for nonlinear response functions provides model-independent re- 
lations between the various components of the RTF thermal n-point functions which can 
be used as consistency checks for approximations. The results derived in this paper can be 
used for elementary or composite bosonic fields with arbitrary interactions. 

For the cases n = 2 and 3 the results from the generalized FDT were shown to reproduce 
known relationships. The new results for the case n = 4 are expected to be useful for 
the derivation of spectral representations for the 4-point function and in a calculation of 
transport coefficients like shear viscosity in scalar field theories. An experimental verification 
of these relations may be possible in applications to the recent developments of picosecond 
pulse techniques and multichannel data acquisition in nonlinear systems p3| . 
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APPENDIX A: RELATION BETWEEN r/a- AND .R/yl-FUNCTIONS 



In this appendix we derive the relation between the r/a vertex functions used in the 



text and the R/A-i unctions introduced by Aurenche and Becherrawy in Ref. I 20fl . As shown 
in 



P0| , the propagator G in the single-time representation can be diagonalized with two 
matrices U and V, 



G = UGV , 



(Al) 



where 



and we choose 



as in Ref. [23 



G 



U 



V 



^ Grr, Gra ^ ^ 



Gar, Gaa 



G re \ 
0, G adv 



Uir, U\a 
U2R, U 2 A 



\ (1 



n 

1, 1 + 71 



V 

f V R1 , V R2 \ (l+n, n ^ 

Vai, Va2 



(A2) 



(A3a) 



(A3b) 



In the CTP formalism the relation between propagator G and the advanced, retarded 
and correlation functions is [DT 



G = Q^GQ 



(A4) 



where 



G 



Q 



g f = q- 1 



Gra, G rr 



^ Qal, Qa2 



0, G 

/ret /~icov 



Qrl, Qr2 
Qla) Qlr ^ 



2r 



V2 
1 

72 



1 / 



\ 



-1, 1 



(A5a) 



(A5b) 



(A5c) 
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Combining Eqs. flA~T|) and ([A4]) leads to 

G = U^Q^GQV- 1 



(A6) 



1. Self energy function 



Substituting the Schwinger- Dyson equation 



G- 1 = G 1 - E 
G- 1 = G 1 - E 



into Eq. ([A6|) we deduce 



where 



V 



^AR, SaA y 

The explicit matrix form of Eq. ( |A8| ) is 



/ W2) F(2) \ 

aa ' ar 

F(2) f( 2 ) , 

V ra ' rr y 



V 



\ 

) 



\ 



(A7a) 
(A7b) 



(At 



(A9) 



^RR, SrA 

£ar, £aa 

rg) + (i + 2n)r( 2 r ), -irf a ) + |(i + 2n)(r( 2 ) - r«) + |(i + 2n) 2 rf r ) 
-2rg», r( 2 )-(i + 2n)r( 2 ,) 

Inserting the FDTs fl62|) and Eq. fl2"0|) into the above expression we get following relations: 



(A10) 



J RR 



, ^RA — ^AR — 0, £_r R — Zj AA 



(All) 



The first two of these equations agree with Eqs. (4.3) and (4.7) in Ref. |20 . 

Please note that and Eaa do not vanish, in contrast to the multi-point cases (c.f. 
Eqs. (|A"ToaD and iM). 
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2. n-point vertex function (n > 3) 



As shown in Ref. [20], the thermal distribution functions which occur in the finite tem- 
perature propagators via the matrices U(k°) and V(k°) can be fully absorbed into the vertex 
functions. To this end one combines the matrices U and V in the propagators (|A1|) with 
the vertex to which G attaches; U(k°) is associated with an outgoing line while V(k°) is 
associated to an incoming line with momentum k. This defines the so-called R/A vertex func- 
tions PD[ . Thus the n-point R/A vertex function with all incoming momenta fci, k-i, ■ ■ ■ , k n 
can be expressed as 

rffijjfci ...**) = v Kiai {kl) ■ ■ ■ vW^)rj£..«.(*i ...*»), (A12) 

where A 1; ...,A n = R, A. f $ An incorporates all temperature dependence of the finite- 
temperature retarded and advanced fields [Q. Inserting Eq. (|24]) into the above equation 
and defining 

P Aiai (k°) = V2V Aiaj (k°)Ql. ai , (A13) 

we obtain the relation between the r /a-functions used in the text and the R/A- functions of 
Aurenche and Becherrawy |20[ : 



f&LiJ*! ■■■*») = 2 1 -P Aiai (A;?) . . . P^M^USh ...kn), (A14) 
where a\, . . . , a n = r, a. The explicit form of Eq. ( |A13| ) reads 



P Ra (k») = 1 , P Rr (k°) = N(k°) = 1 + 2n(k* 
P Aa (k</) = 0, P Ar (k°)=-2. (A15) 

Clearly the R/A vertex functions kn are linear combinations of the r/a vertex functions 
a . We will now show that the R/A vertex functions have simple relations with the 
independent components of the r/a vertex functions after inserting the generalized FDTs. 
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a. 3-point vertex function 



Substituting the FDTs © and Eq. (g) into Eq. (|AT4|) for n=3 we deduce 



f(3) 

1 AAA — 


f (3) _ Q 


(Al6a) 


f(3) 

1 i^AA — 


r(3) 

arr ? 


(A16b) 


f(3) _ 
1 ARA _ 


r(3) 

rar ? 


(A16c) 


f(3) 

1 AAR _ 


F(3) 
rra ) 


(A16d) 


f(3) 

1 ART? — 


~(N 2 + N 3 )V^>, 


(A16e) 


f(3) _ 


-±(N 1 + N 3 )T% r \ 


(A16f) 


f(3) _ 
1 B.RA — 


_I(jv 1 + JV a )rgi*. 


(A16g) 



These relations can be also obtained in the following simpler manner. Noticing 



n 



(*?) 



and substituting Eq. ( [67]) into Eq. ( |A12|) we have 



(P£U) = r£g)(3(3-i( J v 2 + i V3 )rS)r(; 



1\ M /1 N 



(A17) 



(A18) 



Eqs. ( |A16| ) can now be directly read off. It should be noted that the FDTs play an important 
role in getting Eqs. (|A16|) . 

Using f r\ a , T'ara an d ^aar as independent components, Eqs. (|A16|) can be used to 
rewrite Eq. (|67|) as 

" 1 \ / 



f r (3) A = f 
- f 



(3) 
.RAA 



"2 



"3 



, lJ\-(l + n 2 )J\-(l + n 3 )J 
(3) * (1 + n 2 )(l + n 3 ) - n 2 n 3 ( n x 



RAA~ 



'l + ^l) -m 



i + ni ))(-i)(-i, 
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+ f(3) 



ARA I 



ARA~ 



+ f(3) 



n\ \ ( 1 \ / n 3 
l + ni ))\-l) \-(l + n 3 ) / 

(3)* (1 + Wi)(l+W 3 ) -7liW 3 / 1 \ / 

{l + n 2 )-n 2 
( n x \( n 2 \fl" 
AAR {-(l+ ni )) [-(l+r*)) [-1 

(3)* (1 + Tli)(l + n 2 ) -7liW 2 / 1 \ / 

(l + n 3 )-n 3 UlJl 



AAiT 



n 2 



(A19) 



6. J^-point vertex function 
Similar to the 3-point vertex function, using Eqs. (|A12 ) or ( |A14| ) we deduce 



f(4) _ 
1 AAAA — 


f (4) _ Q 
1 RRRR — u J 


(A20a) 


f( 4 ) - 
1 iJiJflA — 


U -l- ni){L -+- ri2jvi + ^3j-p(4)* 


(A20b) 


f(4) 


(l + 7ll )(l + n2 )(l + 7l4 ) 

_ rrar J 
n 3 


(A20c) 


f(4) 

1 RARR — 


(l + 7ii)(l + n 3 )(l + n 4 ) (4 )* 
n 2 


(A20d) 


f(4) 


(l + n2 )(l + n3 )(l + n4 ) 

arrr ' 

ni 


(A20e) 


f(4) 


2 (34) \^ rraa ' iv 3 ± rrra ^ 4-*- rrar,/ ' 


(A20f) 


f(4) 


I/v (13 Yr (4) * -+- A/n"r (4) * + alf^M 

2 (24) ^ rara ' 1 * 2 A rrra ' 2 " 4- 1 - rarr y ) 


(A20g) 


f(4) 


1 /v (14 Vr( 4 )* 4- /v r( 4 )* + n r( 4 M 

2 (23) \ raar ' 1 * 2- 1 - rrar ~r J v 3- 1 - rarr y 5 


(A20h) 


f(4) 

1 AARR — 


Ifr(4) , wrW + /v".r(4) A 

2 y rraa 1 3- 1 - rrra 1 4- 1 - rrar y i 


(A20i) 


f(4) 


Ifr(4) , y\r H 4 ) 4- A^,r( 4 ) 

2 y rara ' rrra ' 4- 1 - rarr y ) 


(A20j) 


f(4) 


-(r (4) + A^ 2 r (4) + Ar,r (4) 1 

2 \ raar 1 ^ rrar 1 J rarr y 5 


(A20k) 


f(4) 

1 AAAi?, — 


_r( 4 ) 

rrra ' 


(A201) 


f(4) 

1 AARA — 


_r( 4 ) 


(A20m) 


f(4) 

1 A.RAA _ 


_r( 4 ) 

rarr ? 


(A20n) 


f(4) 

1 RAAA — 


_r( 4 ) 

arrr 


(A20o) 
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for 4-point vertex function. 

Using f R \ AA , f ARAA , f AARA , F AAAR , ^raaw ^rara> ^rraa as independent compo- 
nents, similar to Eq. ( |A19| ) we can express the 4-point vertex function in the single-time 



representation in terms of column vectors as |32 
' 1 \ / 



r(4) 
1 01020304 



;(4) I 
RAAA \ 



n 2 



n 3 



-1 



[l + n 2 ) \-{l+n 3 ) 



-(4)* (1 + n 2 ) (1 + n 3 ) (1 + n 4 ) - n 2 n 3 n A 

+ 1 RAAA 



+ f(4) 



ARAA I 



m 
l + ni 



1 + m) - rai 



n 3 

l+n 3 ) 



+ r 



( 4 )* (1 +ni)(l + n 3 )(l + n 4 ) -nin 3 n 4 



n 4 \ 
-(l + n 4 V 
/ ni 

V-(i + m 

n 4 \ 

-(l+n 4 V 
/ 1 \ / 



ARAA 



, f (4) I 
+ 1 AAi?A I 



;i+ni: 



(1 + n 2 ) - n 2 
- l + n 2 



-1; 



n 2 

(l + n 2 ) 



■1 V-(l +n 4 ) 



+ r 



(4)* (1 +ni)(l + n 2 )(l + n 4 ) - nin 2 n 4 / 1 



AARA' 



, f (4) I 
+ 1 AAAR 



l + m 



1 + n 3 ) - n 3 
>A- l + n 2 



n 3 



+ r 



(4)* (1 + wi)(l + n 2 )(l + n 3 ) - 72x^2^3 / 1 



AAAR' 



+ r 



(4) 

i?AAi? 



' 1 



(l+n 4 ) 
n 2 \ 

:i+n 2 )J 



n 4 



n 4 



™3 



1 



, f (4)* (1 + w 2 )(l + 7i 3 ) -n 2 ra 3 



+ f (4) 



_RA_RA 



-1 



n 2 



1 



■1, 



+ f(4) 



iJRAA 



m \ 

(1 + ni)(l + n 4 ) - nin 4 + ni) J 

n 4 

-(l+n 4 ), 

m \ 

(1 + m)(i + n 3 ) - nin 3 ^-(1 + m)y 

-(l + n 4 ), 
ni \ 
:i + nx)J 



n 4 



-(4> (1 + n 2 ) (1 + n 4 ) - n 2 n 4 

+ 1 _RA_RA 



™3 



n 4 ) 



' 1 x 



lj(-lj(-(l + n 3 )y 
p(4)* (1 + n 3 )(l + n 4 ) - n 3 ra 4 
(1 + m)(l + n 2 ) - riin 2 



n 2 



(A21) 
51). A 



Please note the much more symmetric structure of this equation compared to Eq. ( 
generalization of Eqs. ( |A19| ) and ( A21|) to n > 5 was recently given in ||32|| . We emphasize 
once more that these representations (in particular the vanishing of R ) rely heavily on 
the validity of the generalized FDTs (or, equivalently, of the KMS condition). Similarly 
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simple expressions can thus not be expected to hold outside of thermal equilibrium. 
From Eqs. ( |A16| ) and ( |A20| ) we can verify the important symmetry relation 



S(n) 
A1...A1 



(fei . . . k n ) = (-1) 



w+ l {i\Aj=R} If ( n )* _ 

II n{kf) Al - A 

{i\A,=A} 



(k\ . . . k n ) , 



(A22) 



where A, = A for A< = R and vice versa. A similar relation was given in Ref. |21 
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